(N 



epi draft 



Branching Ratio and CP Asymmetry of Bg -^ 
i^o*(1430)7r Decays in the PQCD Approach 

Zhi-Qing Zhang 
Department of Physics, Henan University of Technology, 



(N 

o 

fvq ■ Zhengzhou, Henan 450052, P.R.China 



PACS 13 . 25 . Hw - Decays of bottom mesons 
PACS 12 . 38 . Bx - Perturbative calculations 
PACS 14 . 40 . Nd - Bottom mesons 



Abstract - In the two-quark model supposition for 7i'o(1430), the branching ratios and the direct 
r^^ , CP-violating asymmetries for decays B^' — > A'o''(1430)7r°, _K'g'''(1430)7r~ are studied by employing 

CiH' the perturbative QCD factorization approach. We find that although these two decays are both 

tree-dominated, the ratio of their penguin to tree contributions are very different: there is only 
flj . a few percent for the decay B^ — >■ ii'Q'''(1430)7r^, while about 37% in scenario I, even 51% in 

(~| ' scenario II for the decay B^ —>■ Kq^ {14S0)tv'^ . It results that these two decays have very different 

values in the branching ratios and the direct CP asymmetries. The branching ratio of the decay 

B° -> A'*+(1430)7r" is at the order of 10"^ and its direct CP asymmetry is about (20 - 30)% 



'0. 



While for the decay _B^ — > A'o''(1430)7r'', its direct CP-violating asymmetry is very large and about 
90%, but it is difficult to measure it, because the branching ratio for this channel is small and only 
10"'^ order. 
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Introduction. — Along with many scalar mesons found in experiments, more and more 
efforts have been made to study the scalar meson spectrum theoretically [IHZ] • There are two 

k>( ■ typical schemes for their classification [Tl[2] . Scenario I (SI) : the nonet mesons below 1 GeV, 

j_j I including /o(600), /o(980), A'g(800), and ao(980), arc usually viewed as the lowest lying qq 

_C^_- states, while the nonet ones near 1.5 GeV, including /o(1370), /o(1500)//o(1700), A'*(1430), 

and ao(1450), arc suggested as the first excited states. In scenario II (SII), the nonet 
mesons near 1.5 GeV arc treated as qq ground states, while the nonet mesons below 1 
GeV are exotic states beyond the quark model, such as four-quark bound states. In order 
to uncover the inner structures of these scalar mesons, many factorization approaches are 
used to research the B meson decay modes with a final state scalar meson, such as the 
generalized factorization approach [8], QCD factorization approach [Qllllj. and perturbative 
QCD (PQCD) approach [TMS]- 

Whether Kq'^{14:3Q) belongs to the first excited state (scenario I) or the lowest lying 
state (scenario II) is an interesting question, which impacts on its hadronic parameters, 
such as form factor, decay constant and Gegenbauer moment. For example, the form factor 
Fq " [q^) is defined as 

^0 ^'^^ l-a{q^/ml) + b{q^lmir ^' 
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where the parameters a, b are obtained from the fitting procedure among the region < q^ < 
lOGeV^ . They have different values in two scenarios by using the covariant hght-front quark 
model [T7], a = 0.59, b — 0.09 for scenario I and a ~ 0.44, b = 0.05 for scenario II. From the 
potential model calculation, we know that the decay constant and the form factor for sce- 
nario I have opposite signs [lOj . while it is not the case for scenario II. From QCD sum rule 
calculation [18], one can find different masses for the scalar meson sq{q = u, d) under the dif- 
ferent scenario assumptions:m(s5) = 1.410± 0.049 GeV for scenario II and m{sq) > 2.0 GeV 
for scenario I, so the authors considered that scenario II is more favored. B meson decays 
offer a promising opportunity to investigate this question: Here Kq (1430) can be treated as a 
qq state in both scenarios, it is easy to make quantitative predictions in the two-quark model 
supposition, so we would like to use the PQCD approach to calculate the branching ratios 
and the CP-violating asymmetries for decays B° ->■ iV'^°(1430)7r°,is:*+(1430)7r-. Certainly, 
we use the hadronic parameters derived from QCD sum-rule method in our calculations. 
Comparing the theoretical prediction with the future experimental data will indicate infor- 
mation on the structure of A'q°(1430). In the following, T^Tq (1430) is denoted as Kq in some 
places for convenience. 

The perturbative QCD calculation. — Under the two-quark model for the scalar 
meson A'q(1430) supposition, the amplitudes for decays 5° — >• KqTt can be conceptually 
written as the convolution. 



A{Bs -^ K*n) 



dxidx2dx3bidbib2db2b3db3 



•Tr 



C(t)$B. (xi, &i)$if . {x2,b2)<^^ix3,b3)H{x,, h,t)St{x,) e~^^'> ,(2) 



where bi{i = 1,2,3) is the conjugate space coordinate of fc^r, and t is the largest energy 
scale in function H{xi,bi,t). Here Tr denotes the trace over Dirac and color indices, C{t) 
is the Wilson coefficient eval uated a t scale t, which includes the hard dynamics being from 
mw scale down to i ~ ©(a/AMbJ. The function H{ki,k2,k3,t) describes the six-quark 
hard scattering kernel, which consists of the effective four quark operators and a hard gluon 
to connect the spectator quark in the decay. In order to smear the end-point singularity on 
Xi, the jet function St{x) j20) . which comes from the resummation of the double logarithms 
In^Xi, is used. The last term e~^^*'> is the Sudakov form factor which suppresses the soft 
dynamics effectively [H]. So this hard part H can be perturbatively calculated. Here 
Xi{i = 1,2, 3) is the momenta fraction of the antiquark in each meson. There are the same 
conventions with Refs. [15l[T6] in our calculations. 

In the standard model, the related weak effective Hamiltonian H^ff mediating the b ^ d 
type transitions can be written as 



^e// = 



V2 



10 



E ^pi^Vd iCl{^Ji)ol{^x) + C2{ii)oi{p)) - Vt^v;^ J2 ^^(m) ^.(m) 



i=3 



(3) 



where the local four-quark operator Qi{i = 1, ..., 10) and the corresponding Wilson coefficient 



C, can be found in Ref. 
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p(t)b 



Vp(t)d are the CKM matrix elements. 



From the leading order Feynman diagrams for each considered channel, it is easy to get 
the analytic formulas for the amplitudes corresponding to {V — A){V — A), {V — A){V + A) 
and [S-P){S + P) operators, which are similar to those of B ^ /o(980)is:(7r), /o(1500)ii'(7r) 
|151ll6j . We just need to replace some corresponding wave functions. Wilson coefficients, 
and parameters. 

Combining the contributions from different diagrams, the total decay amplitudes for 
these decays can be written as 
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C-i Co 



FeK^ (ae - ^«8) + AfeK* (-^ + ^ 
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PI 
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3G 



X (C5 - ^) + A4k* ^ - MaKS (C3 - -Cg) - FaKS (^4 " ^^lo) 
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+M^k (^5 " J^7) + ^'^^o («4 - ^«io) + F^^, (ae - ^ag) 
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where F^i^jj^* and Afg^^)^* are the n meson emission (annihilation) factorizable contribu- 
tions and nonfactorizable contributions from penguin operators respectively. The upper 
label T denotes the contributions from the tree operators. PI and P2 denote the contribu- 
tions from the {V — A){V + A) and (5 — P){S + P) type operators, respectively. The others 
are the contributions from the {V — A){V — A) type ones. The combinations of the Wilson 
coefficients arc defined as usual: 

ai(^) = C2(m) + ^, a2(M)=Ci(/i) + ^, 

a,(M) = C,(m) + ^^^^, z = 3,5,7,9, 

a,(M) - C,(a.) + ^^M, ^ = 4,6,8,10. (6) 

Numerical results and discussions. — In the two-quark picture, the scalar decay 
constant f^* for the scalar meson A'q can be defined as 



(■^^'ob)l?29i|0) = rriK'jK- 



(7) 



where mx* (j>) is the mass (momentum) of K^ . The light-cone distribution amplitudes for 



the scalar meson Kq can be written as 



{K*o{p)Mz)iq2iO)j\0) 



rm. 



dx ( 



c Jo 



y.{i>^K^,{x) + 7n^.$|.(x) + TUK^ii+i- - l)$^.(x)},z, (8) 



where n+ and n_ are lightlike vectors: n+ = (1,0,0t),?^- = (0,1, Ot). The twist-2 light- 
cone distribution amplitude <^k' can be expanded in the Gegenbauer polynomials: 



*A-*(a;,//) 



/i^g (m) 



6x(l — x) 



i3o(A*)+E^'»(^)^™'(2x-l) 
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As for the twist-3 distribution amplitudes $||-. and $^. , we adopt the asymptotic form: 
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The decay constant fx' and the Gegenbauer moments Bi , B^ of distribution amplitudes for 
A'g(1430) have been calculated in the QCD sum rules [TU], which are listed as 

scenarioI:Si = 0.58 ± 0.07, B3 = -1.2 ± 0.08, /if- = -(300 ± 30)MeV, (11) 

scenarioll :Si = -0.57 ± 0.13, B3 = -0.42 ± 0.22, /k- = (445 ± 50)MeV, (12) 
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Table 1: Decay amplitudes for decays B° -^ A'o''(1430)7r'', /s:*+(1430)7r- (xlO~^GeV^). 







FeK^ 


M,t 


M,KS 


MaKS 


FaK^ 


B'i -^ /Vo*"(1430)7rU (SI) 


-13.4 


-1.2 


-6.4-h4.5i 


0.15-0.08i 


-0.04 -)- 0.08i 


-2.7-6.5i 


B° -^ K*°{U30)Tr° (SII) 


16.5 


1.7 


0.2-h3.9i 


0.01-0.09i 


0.20 + 0. Hi 


0.2 + 8. 6i 


B° -^ if *+ (1430)^- (SI) 


121.3 


0.1 


4.0- 2. 7z 


-0.14 + 0.09i 


0.05 -O.IH 


3.8 + 9. li 


B^ ^ /\*+(1430)7r- (SII) 


-220.7 


0.6 


0.1- 2. 7i 


-0.01-l-O.lli 


-0.27-0.16i 


-0.3-12.2i 



The other input parameters used in the numerical calculations are specified below [23] 



a 

iKfcl 
\Vm\ 



230MeV,MB, = 5.37GeV,Mw = 80.416-6^, 
100° ± 20°, rs, = 1.470 x lO'^^^^ 
3.93 X 10-^Kd = 0.974, 
8.1 X 10~^yt^, = 1.0. 



(13) 
(14) 
(15) 
(16) 



+4.4 



It is noticed that there is dramatic difference for the centra value of CKM angle a between 
the previous and the present PDG values. The present PDG value of a is 89.0 
which is shown in the following figures. 

In the i?s-rest frame, the decay width of B^ 



G|(l-r|,) 



327r77iB, 



/•srQ(1430)7r can be written as 
^\M\^ (17) 



where A4 is the total decay amplitude of each considered decay, which has been given in the 
previous section, and the mass ratio rx* = rriK* /Mb,- M. can be rewritten as 



M = Vu,v:^T - VttV;,P = VubV:, 



1 



ze 



i{-a+S} 



(18) 



where a is the Cabibbo-Kobayashi-Maskawa weak phase angle, and S is the relative strong 
phase between the tree and the penguin amplitudes, which are denoted as "T" and "P," 
respectively. The term z describes the ratio of penguin to tree contributions and is defined 



VtbK 



td 



y^bv:. 



(19) 



If one relates 7W with M., which is the total decay amplitude for the corresponding conju- 
gated decay mode, it is easy to rewrite the decay width F as 



GUI 



'K, 



) 



VubV^aT]^ [l + 2zcosacosJ + z2] 



32TrmBs 
So the CP-averaged branching ratio for each considered decay is defined as 

B^TTBjh, 



(20) 



(21) 



where tb^ is the life time of Bg meson. 

Using the input parameters and the wave functions as specified in this section, it is easy 
to get the values of the factorizable and nonfactorizable amplitudes from the emission and 
annihilation topology diagrams of the considered decays in both scenarios, which are listed 
in Table 1. It is noticed that each penguin amplitude in the table includes the contribu- 
tions from three kinds of operators. The emission factorizable and nonfactorizbale diagram 
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Table 2: CP-averaged branching ratios and direct CP asymmetries of decays _B^ — ^ 

Kq (1430)7r , Kq^ (1430)tv^ . The theoretical errors of branching ratios correspond to the uncer- 
tainties due to variation of (i) the Bs meson shape parameter loi,, (ii) the decay constant /k* , 
(iii) the Gegenbauer moments Bi and -B3 of the scalar meson Kq . While the errors of direct CP 
asymmetries are mainly from (i), (iii) and the CKM angle a = (100 ± 20)°. 



Branching ratios 



(A t; + i.4+U.y+().4N 
1^-3-1.1-0.9-0.4^ 
(A 1+1. 0+1. 0+0. 6>, 
V^--'--0.7-0.1-0.6-' 
/-I 9+0. 5+0. 2+0. In 
V^-^-0.3-0.2-0.1-' 
(o 7+1. 4+0. 9+0. 8\ 
V'^-' -1.0-0.8-0.7-1 



Direct CP asymmetries (in %) 

,+4.8+i. 0+0.0 



i?o ^ A'o*0(1430)7rO (SI) 
B° -> A'*0(1430)7rO (SII) 
B° -^ A'*+(1430)7r- (SI) 
B° -^ A'*+(1430)7r- (SII) 



10-' 
10-7 
10-5 
10-5 



9'J-6_7.3_i.o-8.7 
Qt- r + 1. 2+3. 4+4.0 
^^•"J-8.7-4.1-13.9 
OS 1+4.9+0.7+0.0 
^o-J._4.4-0.7-2.4 
91 n+3-4+2. 1+0.0 
^^■"-3.1-2.5-2.5 



contributions from the tree operators are dominant. In the decay B*^ — >■ ArQ~'"(1430)7r , the 
amplitude F^j^, is enhanced by the large Wilson coefficients C2 + Ci/3, which induce the 
tree operator contributions to be absolutely dominant and the ratio P/T is only 8% in sce- 
nario I, 5.5% in scenario II. Certainly, for the other channel 5° — ?► ArQ"(1430)7r°, the Wilson 
coefficients associated with the amplitude F?^. arc C2/3 -I- Ci, which arc color suppressed. 
We know that the sign of C2/3 is positive while the sign of Ci is negative, which can cancel 
each other mostly. So the amplitude F^j^, is highly suppressed compared with the one for 

the channel S" -^ ArQ~''(1430)7r-. So that the influence from the penguin contributions 
becomes important in the decay 5° — )■ ArQ'^(1430)7r''. The penguin operator contributions 
are mainly from the factorizable emission diagrams and annihilation diagrams, the latter 
are more important and provide a large imagine part, which often induces a large direct 
CP-violating asymmetry (shown in Table 2). From Table 1, one can find that the ratio P/T 
for the decay B'^ — > A'Q°(1430)7r*' is large and about 37% for scenario I, 51% for scenario II. 
Using the amplitudes as specified in Table 1, we can calculate the branching ratios of the 
considered modes, which are listed in Table 2. The uncertainties are mainly from the Bs 
meson shape parameter ujf,, the decay constant /k', the Gegenbauer moments Bi and B3 of 
the scalar meson Kq. From the results, one can find that the branching ratio of the decay 
channel B° -^ A:*+(1430)7r- is about two order lager than that of B° -> A'*°(1430)7r". It is 
because that the former receives a much larger tt" emission factorizable diagram amplitude 
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Fig. 1: The dependence of the branching ratios for B° -^ Ao°(1430)7r'' (a) and B° -^ A*+(1430)7r- 
(b) on the Cabibbo-Kobayashi-Maskawa angle a. The dashed (solid) curves are plotted in scenario 
I (II). The vertical bands show the range of a: 89.0ttt ED- 
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Fig. 2: The dependence of the direct CP- violating asymmetries for B^ -^ ii'o''(1430)7r'' (a) and 
&l — ^ i(rQ + (1430)7r~ (b) on the Cabibbo-Kobayashi-Maskawa angle a. The dashed (solid) curves 
are plotted in scenario I (II). The vertical bands show the range of a: 89.0^42 |24| . 



than the latter. It is interesting to contrast the predictions for decays &l — > A'qtt and those 
for &l — > ii'Q(1430)7r, it seems to exist some similar point: the QCD factorization approach 
predicted the branching ratio of the decay B^ — > KqTt^ was about (1.02J^Qg2) x 10^^ and 
(4.9^35) X 10~^ 22 for the decay 5° — >■ KqTt^. Certainly, the present experimental result 
for the decay Bg — >■ Kq'k~ is about (5.0 ± 1.1) x 10~^ [19|. If the future experimental 
value for the decay B^ — > A'Q'*'(1430)7r^ also falls into 10"^ order, which is less than our 
prediction, it is might because the following two reasons: First, the value of the transition 
from factor _B° — )■ i^Q(1430) predicted by the PQCD approach is larger than the data. 
Second, we only calculate in the leading-order. The higher order contributions might give 
some corrections to the leading order results. Some contributions from the next-to-lcading- 
ordcr (NLO) corrections have been calculated for the decays B — ?> Kn [26], and the results 
show that their branching ratios decrease by about 20% after including the NLO effects. 
These effects might also have an influence on our considered decays. 

The dependence of the branching ratios for the decays 5° — ?► K^^ {1AZ{))ti^ and 5° — > 
/•irQ'''(1430)7r^ on the Cabibbo-Kobayashi-Maskawa angle a is displayed in Fig.l. Compared 
with Fig. 1(a), we know that the branching ratio of the decay B^ — > 7^Q^(1430)7r~ is insen- 
sitive to the angle a, which is shown in Fig. 1(b). From the definition of the CP-averaged 
branching ratio shown in Eq. (|2ip . one can find that if the branching ratio is insensitive to 
the angle a, the coefficient of cosa must be near zero. For the decay B'^ — >■ iirQ~''(1430)7r~, 
the value of z is very small in both scenarios, about 0.1 ~ 0.2, so z^ is smaller and can 
be neglected compared with 1, at the same time, the strong phase angle 5 is about 89° in 
scenario II, so the value of (1 -I- 2z cos a cos (5 -I- z^) is close to 1. Although the strong phase 
angle 5 for scenario I is not so large (about 68.6°), the small value of z makes the branching 
ratio in this scenario also insensitive to a. 

Now, we turn to the evaluations of the direct CP- violating asymmetries of the considered 
decays in the PQCD approach. The direct CP-violating asymmetry can be defined as 



Adir 



\M\^-\M\' 
|7W|2 + |7W|2 



2zsinasin^ 



\ + 2z cos a cos 6 -\- z'^ 



(22) 



Using the calculated ratio z and strong phase 5, it is easy to calculate the numerical values 
of ^^'p (in unit of 10^^) in two scenarios, which are the listed in Table 2. The uncertainties 
are mainly from the Bg meson shape parameter Wf, = 0.5 ± 0.05, the Gegenbauer moments 
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Bi and B3 of the scalar meson Kq, the CKM angle a = (100 ± 20)°. For the decay mode 
iCg"'"(1430)7r~, as discussed above, the value of (1 + 22:cosacos^ + z^) is close to 1 in both 
scenarios, so the corresponding direct CP-violating asymmetry is almost proportional to 
sin a, that is to say it attains the maximum near 90° (shown in Fig. 2b). For the decay mode 
ifg'^(1430)7r°, it receives a very large direct CP-violating asymmetry in both scenarios. It 
is not strange: one can recall that the channel B^ — )■ if "tt^ also receives a large direct CP- 
violating asymmetry, {59 A'^'[2j^)% in the pQCD approach [27], about {4:1.6'^ll'l)% predicted 
by QCDF approach [25] . 

Conclusion. — In this paper, wc calculate the branching ratios and the CP-violating 
asymmetries of decays B° — >• ifQ(1430)7r in the PQCD factorization approach. Using the 
decay constants and light-cone distribution amplitudes derived from QCD sum-rule method, 
we find that although these two decays are both tree-dominated, the ratio of their penguin 
to tree contributions are very different, there is only a few percent for the decay 5° — > 
isrg'''(1430)7r~, while about 37% in scenario I, even 51% in scenario II for the decay B^ — > 
Kq'^ {14:30)tt^ , which results these two decays have very different values in the branching 
ratios and the direct CP asymmetries: The branching ratio of the decay B^ — > i^Tg"'" (1430)7r~ 
is at the order of 10~^, its direct CP asymmetry is about (20 — 30)%. While for the decay 
i3° — > /'iro*'(1430)7r°, its direct CP-violating asymmetry is very large and about 90%, but 
it is difficult to measure it, because the branching ratio for this channel is small and only 
10~^ order. Between these two scenarios, there is less different for the results of the channel 
B° -^ A'*°(1430)7r°, while larger different for those of B° -)■ A:*+(1430)7r-, especially from 
the branching ratio. Because its branching ratio is proportional to the modular square of 
the tree amplitude, which for scenario II is about 3.1 times that of scenario I. So the decay 
B'^ -^ A"Q~''(1430)7r~ can be used to determine the inner structure of A'q (1430) by comparing 
the theoretical prediction with the future experimental data. 
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